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Abstract 

The cosmological particle production in a k = expanding de Sitter uni- 
verse with a Hubble parameter H$ is considered for various values of mass 
or conformal coupling of a free, scalar field. One finds that, for a minimally 
coupled field with mass < m? < 9Hq/4 (except for m 2 = 2Hq), the one- 
mode occupation number grows to unity soon after the physical wavelength 
of the mode becomes larger than the Hubble radius, and afterwards diverges 
as n(t) ~ 0{l){X phys (t) / H Q 1 ) 2v , where v = [9/4 — m 2 /Hq] 1 / 2 . However, 
for a field with m 2 > 9H 2 /A, the occupation number of a mode outside the 
Hubble radius is rapidly oscillating and bounded and does not exceed unity. 
These results, readily generalized for cases of a nonminimal coupling, provide 
a clear argument that the long-wavelength vacuum fluctuations of low-mass 
fields in an inflationary universe do show classical behavior, while those of 
heavy fields do not. The interaction or self-interaction does not appear nec- 
essary for the emergence of classical features, which are entirely due to the 
rapid expansion of the de Sitter background and the upside-down nature of 
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quantum oscillators for modes outside the Hubble radius. 
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I. INTRODUCTION 



Our current understanding of the formation of structures in the Universe rests on the re- 
markable evolution of vacuum fluctuations during the inflationary phase. While stretched by 
the rapid expansion from microscopic to astronomical scales, these fluctuations supposedly 
decohere, and start to behave as (the source of) classical density perturbations. 

To date, several physical pictures and formalisms have been put forward as a means to 
understand this "quantum to classical transition." Most authors approached the problem 
within the general framework of decoherence of a subsystem coupled to the environment. 
Two early applications of this idea to density perturbations in inflationary cosmology are 
those of Sakagami Q, and Brandenberger et al. p|. Unruh and Zurek || examined the 
classical transition in the case of a scalar field coupled to a thermal bath. The general 
criteria for decoherence and the study of cosmological perturbations in the squeezed state 
formalism were given by Albrecht et al. Hu and Calzetta || introduced the technique 
of a closed time path for evaluating the decoherence functional, and studied in depth the 
mutual influence of particle production, dissipation, and decoherence. The same method 
was applied to the problem of density perturbations by Hu et al. ||. Other related works 
are those of Calzetta and Mazzitelli ]7|], who considered decoherence due to particle pro- 
duction, Calzetta and Hu || who studied the decoherence of a mean field by its quantum 
fluctuations, and Calzetta and Gonorazky ||, who investigated fluctuations in a model with 
non-linear coupling. A phenomenological model for the decoherence of density perturbations 
was constructed by Kubotani et al. [fTofl . 

Another category of decoherence studies in inflationary cosmology starts with the pio- 



neering work of Guth and Pi |Tl|, who showed that at late times the wave functions for 
long-wavelength modes behave as classical Gaussian probability distributions. Around the 
same time Lyth |~2| pointed out that the wave packet for massless modes in de Sitter space 
does not spread, and the quantum uncertainty becomes negligible. More recently, Polarski 
and Starobinsky ]13[ and Lesgurges et al. examined particle production in de Sitter 
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space and found a strong squeezing of long wavelength modes which makes their behav- 
ior essentially classical. In all of these cases one looks at the free fields, and the classical 
evolution emerges due to a particular property of the de Sitter background. Polarski and 
Starobinsky [13j call this "decoherence without decoherence," to stress that classical fea- 
tures may emerge even in the absence of an interaction which is deemed essential in usual 
decoherence studies. 

A similar but complementary line of reasoning emerges within the stochastic approach to 
inflation [15| which dramatically exhibits the influence of quantum fluctuations on semiclas- 



sical evolution. Starobinsky |L6| already pointed out the negligible canonical commutator of 
long wavelength modes. Nambu Jl7j constructed a statistical operator for the coarse-grained 



field, and argued that its evolution is effectively classical. In a different picture, Ref. [18 



demonstrated the random walk of a highly peaked wave function for a coarse-grained field, 



in accord with the Langevin equation. Finally, Matacz |19| constructed a new formulation 
of stochastic inflation and investigated the resulting semiclassical evolution. 

All these different physical pictures provide much welcomed complementary insights 
into the problem, and often have vastly different powers in explaining slightly modified or 
more general situations. Some differences in the results and interpretations should not be 
overemphasized at this time. We are still in the phase of collecting new evidence through 
studies of different models, and testing of various decoherence criteria and their comparison. 

In this spirit, I will aprroach here the decoherence of vacuum fluctuations in De Sitter 
space using the most traditional technique of quantum field theory in curved spacetime: par- 



ticle production in an expanding universe [[2(J. One elementarty criteria for the semiclassical 



behavior is to look at the magnitude of the particle occupation number. When this num- 
ber is large, the canonical coordinate and momentum operators approximately commute, 
and the system behaves semiclassicaly. Our goal here is to investigate the decoherence of 
vacuum fluctuations in De Sitter space according to this criteria. Our basic model consists 
of fluctuations of a massive scalar field. We will establish how decoherence depends on the 
mass or conformal coupling of the scalar field. 



Calculations of particle production usually deal with comparing the particle number at 
asymptotically early and late times, or with respect to vacuum states defined in two different 



frames. In the de Sitter case, the calculation of the latter kind has been done by Lapedes [21 



who in this way rederived the existence of a thermal bath seen by an observer in the static 
frame \12\. In contrast, a different quantity is calculated here: the time dependent, one- 
mode occupation number in a spatially flat, expanding frame. This will allows us to follow 
the eventual buildup of (quasi)particle occupation number in the course of the cosmological 
expansion. 

The method to do so is not new, but apparently it has not been applied before to 
the De Sitter case. The probable reason is that since the pioneering works in Ref. p3| most 
efforts concerning fluctuations in inflationary cosmology have been to compute the spectrum 
and the evolution of density perturbations in various models, and to confront these results 
with the observational constraints on the microwave background anisotropy and large scale 
structure. The classical nature of fluctuations is usually assumed, or argued for in passing, 
by invoking Gibbons-Hawking radiation [22], or the horizon in de Sitter space. 

We will return to these and other arguments at the end of this paper. For now, our main 
goal is to consider the simplest of all inflationary scenarios, that of a massive field on a true 
de Sitter background, and to study the conditions under which vacuum fluctuations in this 
model evolve into classical fluctuations. 

The essence of the method is given in Sec. II, and its specification to the de Sitter case 
in Sec. III. Explicit examples for different values of mass are discussed in Sec. IV. Section 
V is the summary of the main results and their implications, including the comparison with 
other works. 



II. PARTICLE PRODUCTION 

Let <p denote a free, minimally coupled scalar field of mass m, and S(rj) the scale factor 
for the Robertson- Walker line element with conformal time. After the substitution = x/S, 
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one can represent the amplitude x(%, v) as the sum or the integral of the modes Xk(v)Qk(%)i 
so that V 2 Q^ = —k 2 Q^. The explicit form of the spatial modes depends on the geometry 
of the three-dimensional space. In either case the field amplitudes are those for a set of 
harmonic oscillators with time dependent frequencies: 

X* + "*fa)Xk = , (2.1) 

with 

u&r,) = m 2 S 2 + k 2 - ^ . (2.2) 

Since, for a free field, the modes decouple and each of them may be considered separately, 
from now on we can drop the subscript k. 

If p denotes the canonical momentum of a mode x we ma Y define the one-mode creation 
and annihilation operators a and a^: 



.= (-J^ x+ ' f , (2.3) 



at =(+)i *i 1+ ' p . (2.4) 
V 2 ^2|w(,)| 

These definitions agree with the standard ones when the frequency is real, and the 
potential is upright. This assures that we will have the same initial vacuum states and 
particle excitations as in the textbook cases. When the frequency is imaginary, and the 
potential upside-down, the presence of \u>\ rather than u> assures that a and a* obey the 
creation/annihilation algebra, and that the particle number still makes sense in that regime. 
One should be aware though that when the frequency is imaginary, the Hamiltonian of the 
excitations is nothing like hou. This is the price that we have to pay for having the sensible 
definition of creation and annihilation operators for an upside-down oscillator, and a sign of 
the very different physics in this case. 

We can now define the one-mode number operator n = a^a. Its rate of change, for our 
case, is found to be 
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^ = ~ H ' (a^ - a 2 ) + l^(at 2 + a 2 ) . (2.5) 
ar? 2 |a;| 2 |a;| 

The two terms in Eq. (|2.5| ) describe the two sources of particle production: the possible 
upside-down shape of the potential, and the possible change in curvature of the potential, 
either upside-down or upright, with time. In particular, if the initial state is the Fock vacuum 
with n = 0, a finite value may develop in time. 

In most of the cases that have been considered before [2(J only the second term con- 



tributes. In some of them the rate of change of frequency is small, and the particle number 
may still be defined as an adiabatic invariant. However, if the frequency is imaginary, there 
will be some particle production due to the first term, and the adiabatic invariance of the 
particle number will not apply. It will be seen in a moment how this works in the de Sitter 
case. 

Although Eq. ( |2.5| ) provides a good illustration of what is going on, it does not offer 
an easy way to calculate the time evolution of the number operator. For this I will use 



the method of Ref. |p4| , with a slight generalization to cover the case of the upside-down 
potential. 

Let Xi an d Xi be the initial conditions at some r]i, and Xiiji) an d X2(ji) are the two 
independent solutions of the oscillator equation (|2.1| ). The general solution may be written 

as 

x(v)=Mv)x t + B(r ] )x f l , (2.6) 
with time dependent coefficients 

m = W^) ' (2 ' 7) 

m = wm ' (2 - 8) 

W( Vi ) = Xi(Vi)X2(Vi)-x'i(Vi)X2(Vi) ■ (2.9) 
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Using A and B one can calculate the off-diagonal Bogoliubov coefficient 612: 



1/2 



1/2 



1 

+ 7;- 



A'{r,) 



2 \uj(rj)u)(r)i)\ 1/2 2 
The one-mode occupation number is then given as 

(0|nfo)|0) = |6 12 (^)| 2 



1/2 



(2.10) 



f2.11^ 



where |0) refers to the Fock vacuum at the initial time r\i. This is the main quantity we are 
interested in. If the occupation number is large, we have that (aa^) = (a^a)[l + 0(l/n)}, 
so, to a high accuracy a and commute, x an d P also approximately commute, and the 
behavior of the mode is essentially classical. This elementary criteria for classicality has 
been used before in the Schroedinger picture by Guth and Pi fll|| . Another Schroedinger 



picture criteria, developed by Lyth [F2]|, will be discussed in the Appendix. 

Note that the canonical commutator has the same form for any choice of creation and 
annihilation operators, so the above criteria for classicality is immune to the usual ambiguity 
of defining particles in curved spacetime. The particle production in this case leads to the 
formation of a classical condensate, and quantum fluctuations are effectively decohered. We 
would like to know what the conditions for this to happen in de Sitter space are so that such 
decohered quantum fluctuations may serve as the source of classical density perturbations 
in inflationary cosmology. 



III. fc = 0DE SITTER EXPANSION 



For spatially flat De Sitter expansion the scale factor evolves as S = —(Horj) -1 , with 
i] G (—oo,0), and Hq = const. Introducing a new independent variable z = —kr), the 
one-mode frequency ( p.2|) becomes 



co 2 (z) = 1 + 



z 2 \m 



(3.1) 



The modes themselves are well known: 



Xi{z) = ^ r V~zH?\z) , X2 {z) = ^^\z) . (3.2) 



#(1,2) 

are Hankel functions, and v = J 9/4 — m 2 /H 2 . These modes are normalized so that 
W(r]) = i at all times. The minus sign in the definition of variable z is to avoid complications 
due to the branch cut of Hankel functions on the negative real axis. Thus, the flow of time 
is from large to small values of z. At early times, when z^ ^> 1, the modes reduce to plane 
waves 

\u{z t )\ = l , Xia ( Z ) = -L e ^— /2- V 4) (33) 

The expression for the off-diagonal Bogoliubov coefficient is now simplified, 

bu(z) = \\^z)\^[B{z) - iA(z)] + ylp[A'(z) + iB'(z)] , (3.4) 

and a straightforward calculation leads to our final expression for the number of particles 
produced in course of the expansion 

"« = ^'^>i 2 + ^Mi 2 -i • < 3 - 5 > 

This expression does not depend on the initial moment Zi, but it correctly reproduces the 
initial vacuum condition n(oo) = in the limit z — > z^ — > oo. 

However, it is the opposite regime, z — > 0, the regime of large expansion, that is of most 
interest, as at that time the wavelength of the mode exceeds the Hubble radius. This is one 
way that inflation solves the horizon problem. 

The behavior of frequency Eq. ( ft.ip in this regime may be seen from Fig. |l|. There 
are two cases: real frequency and bounded classical motion for m 2 > 2Hq, and imaginary 
frequency and "rolling down the hill " for m 2 < 2Hq. The first term in Eq. ( |2.5| ) affects 
the particle number at z < 1 only for m 2 < 2Hq, and is absent for larger masses. 

The moment when these effects turns on is about or after the moment at which the 



fluctuation crosses the Hubble radius. The frequency vanishes at zq = \/2 Jl — m 2 /(2Hq). 
Compared to the moment of crossing the Hubble radius, the two scale factors are related as 



/ m 2 \- 1/2 

S(vo) = S( Vc ) [2 - . (3.6) 

The scalar field which drives inflation is supposed to have a very low mass, m 2 <C H^, 
and in that case, and up to a mass as large as H , the one-mode frequency turns imaginary 
just before the mode crosses the Hubble radius. However, for Hq < m 2 < 2Hq, the "rolling 
down the hill" is postponed to later times, well after the physical wavelength of the mode 
equals the Hubble radius. Either way, it is the evolution at z < 1 that decides the fate of 
the fluctuation. 

The other source of change in the particle number is the variation of frequency with time. 
This is given as 

M' m 2 /H 2 - 2 



(3.7) 



\u\ \z\{m 2 /H$ -2 + z 2 ) 

Figure [2] shows the behavior of this function for several characteristic values of m 2 / Hq. 
The universal growth as z — > is due to the accelerating expansion of the de Sitter universe. 
For all values of mass the rate of change of frequency exceeds unity very soon after the mode 
crosses the Hubble radius. Thus, the second term in Eq. ( ^75|) gives the major contribution 
to the change in particle number at z < 1 for all values of mass. 

There is another regime with a large frequency change in the m 2 < 2Hq case, near point 
z = z where the frequency vanishes. This period lasts just an e-fold of expansion or so, 
before the rate of change in frequency drops close to unity and then starts to grow with 
z — > 0. In this case we have a shrinking size of the phase space cell, rather than its large 
occupancy, so the expected spike in particle number should be considered as a spurious 
effect. 

To proceed, it will be both sufficient and convenient to evaluate Eq. ( |3.5| ) in the approx- 
imation of small z. Consider first the small argument expansion for Hankel functions when 
v is not an integer ||25|| . For the m 2 > cases this means u ^ 0,1. The relevant mode and 



its derivative may be expressed as 

X2 (z) = A v z u+1/2 + B u z 1/2 ~ u + 0(z u+5 / 2 , z~ u+5/2 ) , (3.8) 
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> 2 ( z ) = („ + l/2)A u z»- 1 ' 2 + (-u + 1/2)5, z- 1 " 1 ' 2 . (3.9) 



X 2 



The z/-dependent coefficients are defined as 



1/2 iitv 

The e l7rt/ factors assure that these modes maintain W = i normalization. In two special 
cases the small argument expansions of the modes have logarithmic corrections (which do 
not affect the qualitative behavior of the solutions): 
v = 0, 

X2 {z) = A Q z 1 ' 2 + B,z l ' 2 \n Z - + 0{z*) , (3.11) 

x 2 W = (^ + A)),- 1/2 + Y z ~ 1/2ln l ' (3 ' 12) 

with 

^ 1/2 • C n i , _ 

^ = — + 1^, B^-^ . (3.13) 

v = 1 

X 2(z) = A lZ ^ 2 + B lZ " 2 \n Z - + C.z- 1 ' 2 +0(z 2 ) , (3.14) 

X' 2 (z)= ^-A 1 + Bl )z^ + \b^' 2 \u Z - - ^V 3 / 2 . (3.15) 
tj-1/2 2C — 1 i i 

A = — + Bl = 2^' 671 = ■ (3 - 16) 

C is the Euler constant. 

The behavior of the modes at small z is easy enough to understand. For non-integer 
v there are three different regimes: (a) m 2 / H 2 < 2 or v > 1/2, (b) 2 < m 2 /Hi < 9/4 or 
1/2 > v > 0, (c) m 2 /Hl > 9/4 or i/ G ift. 

In the first case the potential turns upside-down at zq and the field amplitude increases 
without bound, \i ~ \z\ 1 ^ 2 ^ u as z — > 0. In the second case, the potential is always upright 
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but becomes rather shallow at z = 1, below which the small argument expansion starts 
to apply. The z x l 2 ~ v mode dominates again, only now it describes slow settling to the 
minimum. Finally, for v imaginary we have an oscillatory approach to a minimum at a rate 
proportional to z 1//2 . The relative strength of the two modes depends on the mass. In all 
cases the derivative diverges, as the potential is getting steeper and narrower around x = 0. 

This elementary discussion will be sufficient to understand what happens to the particle 
number. Using the preceding small argument expansion, and noting that the frequency may 
be approximated in the same regime as |u;| — > 1 1/4 — v 2 \ l l 2 z~ l , Eq. (|3~5l ) becomes 
1 



n(z) = — 



+ 



1 



1 



v 



V2 r 



\A u \ 2 z v+u * + \B v \ 2 z-( v+v *'> + A v B* v z v ~ v * + A* v B v z v *- v 



v + 1/2| 2 |A,|V +I/ * + I - u + l/2\ 2 \B u \ 2 z 



2|l/4-z/ 2 | 1 /2 
+ {u+ l/2)(-u* + l/2)A' v Blz v - 



+ {-v + l/2){v* + 1/2)^^*-"] . 



(3.17) 



In the two special cases when m 2 > and v is an integer, we have, at z <C 1, 



n(z) 



1 1 



if i/ = 1 



(3.18) 



n[z) 



2tt 



In- 



if v = 



(3.19) 



We can now see how the mass of the scalar field affects the amount of particle production. 



IV. THE EFFECT OF MASS 



Since from very early on it has been realized that the inflaton field needs to be very 



weekly (self-) coupled, or nearly massless |?3] most of the work on fluctuations in inflationary 
cosmology concentrated on those cases. There are very few results about the behavior of the 
fields with mass of the order of the Hubble parameter, and almost none for heavier fields. 
Yet, there have been clear indications that the behavior of massive scalars in De Sitter space 
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very much depends on the mass. Equation ( ft.lTQ brings this expectation to the full force, 
and immediately yields our main result: The particle number for heavy fields in De Sitter 
space, with m 2 / Hq > 9/4, is an oscillatory function at times later than the Hubble radius 
crossing time. There is some particle production, but there are no classical condensates 
with n>l. In contrast, the particle number for medium and light fields in De Sitter space, 
with < m 2 / Hq < 9/4, except m 2 = 2Hq, diverges at times later than the Hubble radius 
crossing time. There is a significant particle production, leading to the formation of classical 
condensates with n 1. 

The m 2 = 2Hq case is conformally related to flat spacetime. The frequency vanishes 
for this choice of mass, the modes are plane waves at all times, and there is no particle 



production. This is one version of a classic result by Parker |20 

Let us now illustrate these results by taking a closer look at the specific cases. 

(a) m 2 1 Hq = —4, v = 5/2. This kind of model was analyzed by Guth and Pi [Tl| 



as it has direct bearing for the generation of density perturbations in the new inflationary 
scenario. The \m 2 \/HQ ratio is supposed to be very small, but for illustrative purposes I will 
take it here to be large, so that the parameter v is a half-integer, and the modes are given 
through elementary functions. We can then calculate the particle number at all times. This 
is not a limitation, since, according to Eq. Q3.17D , the behavior of the particle number is 
qualitatively similar for all v > 0. We find 



X2iz) = 7=2 



3 . / 3 

z V z z 



u 2 (z) = 1 - A (4.1) 



z 2 



and 



,s 2-3/z 2 + 36/z 6 1 |2 c , ll2 (\ 3 9\ 1 „ n , 

The behavior of this function is shown in Fig. [| The oscillator turns upside down 
while the mode is still within the Hubble radius, and soon after the particle number rapidly 
increases, approaching quickly a z~ 2v growth for wavelength longer than the Hubble radius. 

(b) m 2 = 0, v = 3/2. This model is the starting point for most inflationary scenarios. 
The modes are again given through elementary functions 
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X2[Z) 

and the particle number Eq. 
The result is 

n(z) = ( 



u 2 (z) = 1 - - 



(4.3) 



may be computed through all times, not just at small z. 



)- + -K 
2 4 1 

+ 



1 ! f--l 



1 1 

- + - 

z z 6 



(4.4) 



4|z 2 -2|V2 

The divergence at Zq = ^/2, is due to the vanishing frequency at this point. It is 
represented by the spike on Fig. [| n quickly returns to value below unity as soon as the 
mode crosses the Hubble radius, and then the rapid particle production starts in earnest. For 
long- wavelength modes, well outside the Hubble radius (z <C 1), we have that the occupation 
number diverges as z~ 3 . This may be written as 



n{rj) 



(4.5) 



4v/2V H o' 

(c) < m 2 /Hl < 2, v G (1/2,3/2]. The result for v = 1 is given before Eq. (|3~T8|) . For 
v 7^ 1 the expression for particle number at small z becomes 

~u - 1/2" 



n(z) = v 



1 V 2 2 2i/ - 2 



7T 



Y{v)fz 



2 „-2v 



(4.6) 



_z/ + l/2_ 

Both the growth rate and the overall normalization of this function decrease with the 
mass square increasing from to a zero production value of 2Hq. The behavior smoothly 
changes as v decreases, and it is well represented by the v = 1 case in Fig. |3|. 

(d) 2 < m 2 / Hq < 9/4, v e [0, 1/2). Unlike the preceding case, the frequency now stays 
real and the potential is upright at all times, including times later than the Hubble crossing 
time. Nevertheless, as Eq. ( |3.17| ) shows, the particle number diverges as z — > 0: 

"1/2 - v 



n[z) 



V 2 2 2i/ -3 



-2v 



7T 



(4.7) 



1/2 + p_ 

One example of this kind, with m 2 / Hq = 35/16 and v = 1/4, is illustrated on Fig. |3|. 
The behavior varies smoothly with the increasing mass, until v = 0, given separately in 
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( p,19| ), when the growth drops from a power law to the logarithm. Since the frequency is 
real, the particle production in all of these cases must be entirely due to the rapid change 
of frequency at small z. The relative shallowness of the potential early on (just after the 
Hubble radius crossing) probably plays some role as well. 

(e) m 2 1 Hq > 9/4, v G ilZ. In this case the behavior of n(z) at z <C 1 is purely oscillatory. 
An explicit substitution of v = to Eq. ( [3.17|) gives the following result: 



n(z) = A + A c cos[2|i/|ln(>/2) - 2$ r ] + A s sin [2\v\ ln(z/2) - 2$ r ] . (4.8) 



$r is the phase of r(i|i/|) [pq| . The three coefficients are found to be 



A = (-)- + £ , ' coth(|z/|vr) , (4.9) 





1 








2 


3/4- 


z/| 2 


2\v\ 



Ac = -^UT-^ TTTTT^TT^ e"^ , (4.10) 



1 + COth(|l/|7T) M 

s_l J 4(1/4+ |^| 2 )V2 • ^ AA J 

The period of the oscillations is zi/^ 2 = exp[7r/|z/|]. Figure || shows two such cases. As 
in the case (d), the potential is upright, and the only source of particle production is the 
rapid variation in frequency. However, unlike the previous case, the steeper potential does 
not allow for a large amplitude in x, and the change in frequency does not lead to a large 
particle number, merely to oscillations in its value. The particle number is bounded by, 
n(z) < f = A + (A C + A s ) l/2 . If we denote x = \u\, this amplitude may be expressed as 

1 (1/4 + X 2 ) 1 ' 2 , , . Rix) 1 ' 2 . . . . 
f(x =-» + - o cothTrx + .y ,XG l,oo , 4.12 

2 2x 8 smh(7rxj 

where 

R(x) = * 4 + 5x 2 /2 + 9/l6 
[) ~ x 2 (x 2 + l/A) ' { > 

As Fig. |4] shows, this amplitude is a remarkably smooth function of \v\. The growth as 
^ — > is due to the inaccuracy of the employed small 2 expansion for noninteger v near 
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integer value v = 0. The transition from the bounded and oscillatory behavior of the mode 
at v G iR, to unbounded logarithmic growth when v — cannot be described accurately 
by the lowest order expansion in z, and it is simulated as the divergence of T(z/) as v — > 0. 
However, for |z/| > 1 the expansion may be trusted and, as Fig. [| shows, the particle number 
is indeed bounded for all values of mass. The amplitude of the oscillations tends to zero as 
Hq/itl 2 in the limit of very high mass. 

(f) nonminimal coupling. The effect of adding an £i?0 2 /2 term is particularly simple in 
case of a k = de Sitter background: in all the expressions above one should replace m 2 with 
m 2 + 12^Hq. The behavior of the occupation number as parametrized by different values of 
v is the same as before, except that now v 2 = 9/4 — m 2 / Hq — 12£, and classical condensates 
form whenever this takes positive values or zero. Figure |5] shows the characteristic domains. 
The basic physics of the phenomena is unchanged. 

V. CONCLUSIONS 

As we have seen, the ratio of the two scales, the Hubble radius Hq 1 and the Compton 
wavelength m -1 , determines the onset of two phenomena in de Sitter space: the rapid rate of 
change in frequency of quantum oscillators and/or their turning into the upside-down shape. 
Together or separately these two phenomena lead to the formation of classical condensates 
and effectively classical behavior of quantum fluctuations. 

In this picture parameter Hq 1 is important not because of its causal property, but because 
it is the characteristic scale for spacetime curvature. The decoherence scale varies from 
roughly Hq 1 in the case of a massless field to infinity as mass of the field increases towards 
m 2 = 2Hq. The rapid growth in particle number is due to both the upside-down shape of 
the quantum oscillators and the rapid change in their frequency. For masses in the narrow 
range m 2 G (2Hq,9Hq/4) the decoherence scale is again just a bit larger than the Hubble 
radius. The potential is upright in this case, and the formation of the condensate is entirely 
due to rapid expansion. Finally, for fluctuations with even larger mass the particle number 
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is always bounded. These fluctuations do not decohere, and cannot serve as the source of 
perturbations that seed galaxies. 

The physical picture used here is just another side of the well known phenomena of 
parametric amplification and squeezing of quantum fluctuations |27|, that played an 
important role early on in the studies of particle production and density perturbations. 
Two new moments covered here are the attention to decoherence through the formation 
of classical condensates, and the role of mass in the efficiency of particle production and 
decoherence. 

It is a somewhat curious result that the n(z) ~ z~ 2u growth of the occupation number 
for low and intermediate masses is the same as the behavior of the dispersion ratio used by 



Guth and Pi [|TTJ to characterize the late time classical behavior in similar models. However, 
since they carried out their calculations for a model with a potential of Coleman- Weinberg 
shape, the oscillators were upside-down because of the curvature of that potential, and it 
was not clear from Ref. JTTJ] whether the classical behavior would emerge if the potential in 
the field Lagrangian was upright. Besides elucidating this role of mass, we obtained here a 



somewhat unexpected addendum to Ref. the rate of change of frequency may as well 



be just as important for decoherence of low mass modes as the turning of the oscillators to 
an upside-down shape. In fact, for intermediate mass fields, case (d) above, this is the only 
source of decoherence. 

Another curiosity has to do with the kinematic argument that rapid expansion leads 
to the accumulation of modes with low, and almost indistinguishable frequencies. This 
"infrared divergence" supposedly leads to a large occupation number for long- wavelength 
modes, hence to their highly classical behavior. This is a puzzling argument, since it does not 
depend on the mass of the field, so one would expect the formation of a classical condensate 
in all cases. The answer is that while the occupation number in the massless case indeed 
looks as if it can be understood by a kinematic argument n ~ (\ p h ys / Hq 1 ) 3 for a finite mass 
the exponent is smaller, and precludes kinematic interpretation. 

The decoherence argument which is often considered to be the most trustworthy is the 
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demonstration of the effectively diagonal form of the reduced density matrix obtained by 
tracing out some uninteresting degree (s) of freedom other than inflaton field. This general 
scheme appears to be the most satisfactory approach to the classic "measurement problem" 
in quantum mechanics, and its application to the creation of density perturbations in infla- 
tionary cosmology has been investigated by several authors, e.g., Brandenberger et al. 0, 
Sakagami fl|], and Hu et al. ||. A common conclusion to all these approaches is that the 
decoherence is proportional to the strength of the coupling, and it is supposed to disappear 
if the coupling goes to zero. The analysis of free models done here shows that in the de 
Sitter background there will be some decoherence which does not depend on the magnitude 
of the coupling constant. The interaction or self-interaction may be sufficient but not a 
necessary condition for the decoherence. A similar conclusion may be reached by utilizing 
the stochastic approach of Starobinsky |L5[] for the analysis of long-wavelength quantum 
fluctuations in an inflationary universe [ljj , . 



APPENDIX: COMPARISON WITH LYTH'S CRITERIA FOR DECOHERENCE 

A plausible physical argument for the decoherence of quantum fluctuations in an infla- 
tionary universe was pointed out by Lyth. He considered a massless, minimally coupled 



scalar field, initially in the state of minimal uncertainty. Subsequently, the wave function 
of this time-dependent oscillator will spread. By treating this spread approximately as in 
the case of a free particle in flat spacetime, Lyth was able to show that times taken for the 
width to reach values that are large compared to the initial spread are themselves quite large 
compared to the Hubble time. Due to technical restrictions the argument applies only for a 
period of a few Hubble times after the fluctuation crosses the Hubble radius, but, as Lyth 
remarked, "it seems plausible [...] that if the wave packet is already extremely classical when 
this stage is reached, it will remain so for 60 or so Hubble times that the field survives." 

Indeed, this argument is quite robust, and for a massless case it is equivalent to the 
criteria used here. To see this, recall that Lyth uses the field amplitude before rescalling 
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4> = X2/S. In the massless case, xi is given by Eq. (fO|), so we have 

0(77) = ^Ai - krj) e~ ikn ■ (Al) 

Up to a constant, this mode agrees with Eq. (53) in Ref. [|1^ . 

At late times, \2 diverges, while freezes. The latter property has often been used 
in early studies as a heuristic argument for the scale invariant nature of fluctuations in 
an inflationary universe. The expectation value for the Heisenberg operator 2 will also 
approach a constant. In the Schroedinger picture this property translates into the freezing 
out of the width for the one-mode wave function. Therefore, the "rolling down the hill" in 
the x representation, which is responsible for the growth of the particle number, is equivalent 
in the massless case to the Lyth's criteria (in (f) representation) of a wave packet which does 
not spread appreciably. The two criteria for classical behavior are in agreement. 
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FIGURES 



FIG. 1. The evolution of the one- mode frequency squared in de Sitter space, Eq. (3.1), for 
different values of mass parameter a = m 2 /Hq. In all cases the frequency starts to change rapidly 
after the physical wavelength of the mode exceeds the Hubble radius. In the same regime the 
one- mode quantum oscillators turn to an upside-down shape for a < 2. 

FIG. 2. The evolution of a relative change of the one-mode frequency in de Sitter space, Eq. 
( |3.1| ), for different values of mass. For a < 2 the frequency passes through zero, and there is a 
spurious singularity in |u//u;|. For all values of mass the relative change in frequency exceeds unity 
shortly after the wavelength surpasses the Hubble radius. 

FIG. 3. The evolution of a one-mode occupation number in De Sitter space for different values 
of mass, the main result of this work. The "Coleman- Weinberg " type case with a large mass, 
\a\ = 4, Eq. ( |4.2| ), and the fiducial m? = case, Eq. ( |4.4| ), are plotted for the whole range, 
others only for \kr)\ <C 1, when the mode is well outside the Hubble radius. The divergence at 
kr] = -2 1 / 2 for a = is a spurious effect, due to to 2 passing through zero as the oscillator turns 
into an upside-down shape. The corresponding divergence in the a = —4 case is just outside the 
range of the figure. The a = 5/4 case corresponds to v — 1, and, as a special case, it is given by 
Eq. ( HUp . In a = 35/16, or v = 1/4 case, the particle number may be evaluated from Eq. ( p~7|) 



as n(z) ~ 0.33/z 1//2 . Unlike the preceding two cases the quantum oscillator has upside-right shape 
at \krj\ <C 1, but the particle number nevertheless diverges due to the rapid change in frequency. 
The two cases with imaginary v are a = 3 or v = z3 1 / 2 /2, and a = 25/4 or v = 2i. Both are given 
through Eq. ( |4.8| ). The particle number is bounded, oscillatory function of In \ krj\, and there is no 
formation of classical condensates. 
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FIG. 4. The maximal amplitude of the particle number for heavy fields, m 2 > 9Hq/A, Eq. 
( 4.12] ), as function of mass, x = \v\. This surprisingly smooth function shows that in course of 
its oscillations the particle number never exceeds unity, so the vacuum fluctuations of heavy fields 
in de Sitter space do not have a classical behavior even for wavelengths that exceed the Hubble 
radius. 

FIG. 5. Formation of classical condensates and the decoherence of vacuum fluctuations in De 
Sitter space as function of both the mass and the nonminimal coupling. Condensates form not just 
for lo 2 < 1 and v E TZ, but also for a narrow set of values between u = and uj 2 = 1 (or v = 1/2). 
See part f, Section 4. 
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